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We present a systematic treatment of the initial conditions and evolution of cosmological per- 
turbations in a universe containing photons, baryons, neutrinos, cold dark matter, and a scalar 
quintessence field. By formulating the evolution in terms of a differential equation involving 
a matrix acting on a vector comprised of the perturbation variables, we can use the familiar 
language of eigenvalues and eigenvectors. As the largest eigenvalue of the evolution matrix is 
fourfold degenerate, it follows that there are four dominant modes with non-diverging gravi- 
tational potential at early times, corresponding to adiabatic, cold dark matter isocurvature, 
baryon isocurvature and neutrino isocurvature perturbations. We conclude that quintessence 
does not lead to an additional independent mode. 



Introduction 

We present here a systematic treatment of initial conditions for quintessence models. The 
universe we will consider contains photons, baryons, neutrinos, cold dark matter and a scalar 
quintessence field. We will formulate the evolution equations for the perturbation variables 



U contains all perturbation variables and the matrix A(x) encodes the evolution equations. 
In doing so, we relate the problem of finding initial conditions and dominant modes to the 
familiar language of eigenvalues and eigenvectors. This formulation makes "mode-accounting" 
transparent by counting the degeneracy of the largest eigenvalue. We find four dominant modes 
that remain regular at early times. For physical reasons, we choose a basis using adiabatic, 
CDM isocurvature, baryon isocurvature and neutrino isocurvature initial conditions. 

The Different Modes 

In the following we adopt the gauge-invariant approach as devised by Bardeen^. For a more 
detailed derivation of the perturbation equations see Doran et alH It turns out that the evo- 
lution is best described as a function of x = kr, where r is the conformal time and k the 
comoving wavenumber of the mode. We assume that at early times, the universe expands as 
if radiation dominated. Assuming tracking quintessence we obtain the following set of equations : 



as a first order differential matrix equation: 
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U = A(x)U, (1) 



where the vector 
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with the gauge-invariant Newtonian potential We denote the derivative d/dlnx with a 
prime. The gauge-invariant ene rgy de nsity contrasts A a , the velocities V a and the shear 11^ are 
the ones found in the literature^ ' ' \ except that we factor out powers of x from the velocity 
and shear defining V = V/x and Ilj, = x~ 2 U u . The index a runs over the five species in our 
equations, quintessence is assigned the subscript q. We assume tight coupling between photons 
and baryons. The equation of state w = p/p takes on the values w c = Wb = 0, u) 7 = w u = 1/3 
and w q is left as a free parameter. 

We define the perturbation vector as jj t = (A c , V c , A 7 , V 7 , A b , A v , V u , U u , A q , V q ). (13) 

The matrix A{x) can easily be read off from equations ip])-(|llj). This enables us to discuss the 
problem of specifying initial conditions in a systematic way. 

The initial conditions are specified for modes well outside the horizon, i.e. i < 1, In this 
case, the r.h.s. of equations (j2J, ©, © and ® can be neglected, provided V a does not diverge 
oc x~ 2 or faster for x 2 — > 0. 

The general solution to Equation Q in the (ideal) case of a truly constant A would be 

^) = E c <(£) Alc/W ' ( 14 ) 

where U^' are the eigenvectors of A with eigenvalue Aj and the time independent coefficients c\ 
specify the initial contribution of 17 ^ towards a general perturbation U. As time progresses, 
components corresponding to the largest eigenvalues A, will dominate. Compared to these "dom- 
inant" modes, initial contributions in the direction of eigenvectors with smaller Re(Aj) decay. 
In our case, the characteristic polynomial of A(x) indeed has a fourfold degenerate eigenvalue 
A = in the limit x 2 — > 0, the other six remaining eigenvalues are negative. We therefore need 
to solve A(x)TJ = which is equivalent to setting the l.h.s. of Equations (|2|)- (|11|) equal to zero 
and using Q c = = x 2 = 0. Then Equations @, 0), © and Q are automatically satis- 
fied (provided V a does not diverge oc x~ 2 or faster), and Equations (|3|).([5|).(|8|)- (jllj) yield non- 
trivial constraints for the components of U : 
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1/4A 7 - % + fi„n„ + 2^ = (16) 



1/4A„ - V v + n u fl u + 2* = (17) 8/5K - 2lX = (18) 



3J2 v n i/ +A,/(l+t(;-)+3Fg+3* = (19) 3fi 1 ,n i , + A g /(l + u>,) + F g + 4* = (20) 



Following the existing literature, we use the gauge-invariant entropy perturbation 1 ^ between 
two species a and j3, as well as the gauge-invariant curvature perturbatio n qn hyp er-surfaces of 
uniform energy density of species a in order to classify the physical modes 



Sa-p = T^-T^, (21) Ca=(H L + ±H T ) + , 6pa , . (22) 

In our variables, these expressions take on the manifestly gauge- invariant form 

^ 3(1 + Wa ) ' ^ E a 3(i+™ Q )o; 1 ; 

The first (rather intuitive) perturbations one would try to find are adiabatic perturbations, 
which are specified by the adiabaticity conditions S a -p = for all pairs of components, i.e. 

A„ = A 7 = ^A c = ^A b , (24) 

Using the six constraint Equations (|15 )l -(|20 )l . we obtain the adiabatic mode. (Due to limited 
leneht of this article we cannot quote the full results, we therefore refer the reader to Doran et. 
al. ™for more details.) We also conclude that quintessence is automatically adiabatic if CDM, 
baryons, neutrinos and radiation are adiabatic, independent of the quintessence model for as 
long as we are in the tracking regime. 

Let us next consider the neutrino isocurvature mode. For this, we require that CDM, baryons 
and radiation are adiabatic, while S„.~ ^ and that the gauge-invariant curvature perturbation 
vanishes: 

C = 0, A c = A b = ^A 7 . (25) 

Using this and Equations p5 |) -l|2U |) leads to the neutrino-isocurvature perturbation. 

It is important to note that we did not require quintessence to be adiabatic. One can see 
from the neutrino isocurvature vector that A q = 0, and as a consequence quintessence is not 
adiabatic with respect to either neutrinos, radiation, baryons or CDM. Hence, we could just as 
well have labeled this vector "quintessence isocurvature". 

The CDM isocurvature mode is characterized by S c .~ ^ 0, Q = and adiabaticity between 
photons, neutrinos and baryons. Similarly, for the baryon isocurvature mode we require Sb-y ^ 0, 

C = o. 

The adiabatic, CDM isocurvature, baryon isocurvature and neutrino isocurvature- vector 
are linearly independent. We have therefore identified four modes corresponding to the fourfold 
degenerate eigenvalue zero of A(x). These four vectors span the subspace of dominant modes in 
the super-horizon limit. Arbitrary initial perturbations may therefore be represented by project- 
ing a perturbation vector U at initial time into the subspace spanned by the four aforementioned 
vectors, as this is the part of the initial perturbations which will dominate as time progresses. 

We use a modified version of CMBEASY ^^^ to compute CMB spectra corresponding to 
different initial conditions for an early quintessence cosmology with parameters as in model A 
We set the spectral index of the isocurvature modes identical to the spectral index of the 
pure adiabatic mode, n s = 0.99. Comparison with the WMAP data in the same figure shows 
that non-adiabatic initial perturbations are strongly constrained. 
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Figure 1: CMB Temperature spectra as a function of multipole / in an early quintessence cosmology. The pure 
adiabatic (AD), CDM isocurvature (CI), neutrino isocurvature (NI) mode and three different combinations of these 
dominant modes are plotted. For comparison with experimental data we also give the WMAP measurements of 
the CMB^l The spectrum of the pure baryon isocurvature mode is essentially identical to that of the pure CDM 
isocurvature mode. All spectra have been normalized to the same power at / = 10 

Conclusion 

We have investigated perturbations in a radiation-dominated universe containing quintessence, 
CDM, neutrinos, radiation and baryons in the tight coupling limit. The perturbation evolution 
has been expressed as a differential equation involving a matrix acting on a vector comprised 
of the perturbation variables. This formulation leads to a systematic determination of the ini- 
tial conditions. In particular, we find that due to the presence of tracking scalar quintessence 
no additional dominant mode is introduced. This fact is beautifully transparent in the matrix 
language. In total, we find four dominant modes and choose them as adiabatic, CDM isocur- 
vature, baryon isocurvature and neutrino isocurvature. For the neutrino isocurvature mode, 
quintessence automatically is forced to non-adiabaticity. Hence, we could have as well labeled 
the neutrino isocurvature mode as quintessence isocurvature. To demonstrate the influence on 
the cosmic microwave background anisotropy spectrum, we have calculated spectra for all modes. 
A detailed study may provide ways to put additional constraints on quintessence models or tell 
us more about the initial perturbations after inflation. 
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